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INVLSTItiATION  OF  DISCRLTL  FUNCTION  TFCHNOLOOY 
I OR  TOPOGRAPHIC  SCIFNCtS 

INTRODUCTION  ■ The  purpose  of  this  Work  Unit  was  to  determine  the 
Feasibility  and  potential  merits  of  discrete  function  technology  for  applications  in 
Objectives  the  topographic  sciences,  emphasizing  the  extraction  of  carto- 

graphic features  from  images.  Specific  objectives  of  this  study 
were  to 


1.  Research  the  use  of  various  orthogonal  sets  of  functions  for 
decomposing  topographic  images  in  spectral  components. 

2.  Determine  whether  such  decompositions  or  transforms  are 
useful  for  the  extraction  of  cartographic  anil  terrain  features. 
2.  Compare  the  use  of  Fourier.  Block,  Walsh,  and  Haar  trans- 
forms of  images  for  machine  feature  extraction. 

4.  Research  the  potential  of  devices  that  may  facilitate  direct 
discrete  transforms  using  electro-optical  technology. 

5.  Research  other  applications  of  discrete  functions  to  the 
topographic  sciences. 


Image  analysis  for  the  purpose  of  enhancement,  restoration,  coding,  and  data  and 
feature  extraction  emerged  as  a new  field  of  research  based  primarily  on  the  newly 
Background  developed  technologies  for  television  CI  V)  and  digital  com- 

puters. 

In  TV,  the  image  of  the  TV  camera  is  scanned  and  converted  into  a time- 
dependent  electrical  signal.  This  signal  is  transmitted  through  communication 
channels  to  the  receiver  and  restored  to  the  image.  To  determine  optimum  channel 
bandwidth,  filter  performance,  and  possible  bandwidth  compression,  the  signal 
would  be  transformed  from  the  time  domain  into  the  frequency  domain.  For 
the  most  part,  this  transformation  is  accomplished  in  communeiations  by  the 
application  of  the  Fourier  transform  and  is  also  called  the  decomposition  of  the 
signal  in  its  spectral  components.  If  the  spectral  components  of  a signal  are  known, 
suitable  fillers  may  be  designed  that  reject  most  other  signals  and  noise  and  pass 
only  the  signal  lor  which  the  filter  is  designed.  Thus,  filters  may  be  used  for  sig- 
nal processing  and  extraction.  A IV  signal  consists  of  a carrier,  video  modulation 
representing  the  image  content,  sound  modulation,  and  synchronization  signals. 
The  various  components  of  the  TV  signal  are  separated,  extracted,  and  processed 
by  filters  having  special  purpose  spectral  response  characteristics.  Certain  pro- 
cessing functions  of  a TV  set,  namely  brightness,  contrast,  hue,  tint,  may  be 
adjusted  by  the  viewer  to  his  satisfaction.  These  processing  functions  are  also  used 
in  image  analysis. 

From  digital  computer  technology  emerged  the  simulation  of  signal  processing 
capabilities  and  systems.  The  large  storage  capacity  of  digital  computers  and  the 
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legibility  of  digitizing  images  in  a way  that  they  ean  he  stored  and  processed 
in  a computer  led  to  the  techniques  of  digital  image  processing  and  visual  pattern 
recognition. 

Research  was  conducted  on  the  discrete  lunction  technology  as  used  in  the  analy- 
sis and  leature  extraction  ol  topographic  images,  in  representing  I he  earth's  gravity 
field,  and  in  remote  sensing  applications  ol  discrete  electro- 
scope and  Approach  l,,aSnctic  wave  lorms.  I he  theory,  generation,  and  properties 
ol  potentially  uselul  discrete  functions  are  explained  by  ex- 
amples. I he  applications  to  topographic  sciences  are  investi- 
gated and  also  demonstrated  by  selected  examples.  In  addition,  various  hardware 
approaches  lor  implementation  ol  discrete  lunction  technology  are  critically 
examined,  and  some  concepts  ol  machine  leature  extraction  from  topographic 
images  are  presented. 


THEORY  AND  PROPERTIES 
OF  DISCRETE  F UNC  TIONS 


Basic  Concepts 


■ A rule  of  correspondence  that  associates  a 

real  or  complex  number,  y = f(xj,  xs xn). 

with  each  given  set  of  real  or  complex 
numbers  Xj,  xs,.  . x„  is  called  a lunction 
ol  x | , Xi,  . . . xn.  The  lunction  can  be  expressed  analytically  or  by  an  n-dimen- 
sional  table.  I he  arguments  may  assume  all  real  and  complex  numbers  between 
-°°  and  +«>  and/or  portions  thereof.  Many  of  the  functions  commonly  treated  in 
applied  mathematics  are  continuous  and  have  finite  derivatives. 


I he  y-values  ol  discrete  lunction  are  frequently  real  and  vary  only  in  discrete 
steps.  Many  discrete  functions  of  interest  assume  only  the  values 

- 1 , 0,  + 1 or  just  0 and  I . 


II  the  discrete  function  assumes  only  two  values,  such  as  0 and  I or  -I  and  +1, 
the  lunction  is  also  called  a binary  function.  Most  discrete  functions  have  in- 
teger values  and  rarely  exceed  1000.  Discrete  functions  arc  not  necessarily  con- 
tinuous and  possess  no  derivatives  in  the  normal  sense. 

Discrete  functions  used  in  computers  and  communication  equipment  include 
the  output  ol  analog-to-digital  converters,  digital  voltmeters,  digital  clocks, 
maximum  length  code  generators,  and  teletype  machines.  This  report  deals  with 
discrete  (unctions  as  they  pertain  to  the  time,  the  spatial  coordinates  x and  y 
of  a plane,  or  the  three-dimensional  space  coordinates. 

Signals  for  TV  have  two  space  variables  and  one  time  variable.  The  space  variables 
are  derived  from  the  image  of  the  TV  camera  by  a scanner  in  the  case  of  digital 
Block  Pulse  TV.  The  scanner  operation  may  be  described  by  a two-dimen- 

Function  sional  block  pulse  function. 
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f igure  1 (.rapine  Representation  of  the 
Block  Pulse  1 unctions 
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i I wo-Oimensional  Block  Boise  I*  unctions 


Bln  (k,x)  bln  (t,y) 

K 0,1, 2,3,4, 5,6, 7, 
v 0,1, 2, .1,4, 1,6, 7, 


b o (A,  k)  r,  0 / 

n * -• 
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The  block  pulse  function  blofk.x)  is 
defined  in  the  interval  -Vi  ( x < +Vi.  Ii 
has  the  values  0 and  +1 . I he  pulse  shilts 
from  left  to  right  within  the  interval. 

The  parameter  K = 0,  1,2, ,n  def  ines 

the  number  of  shif  ts  in  equal  steps  with- 
in the  interval.  figure  I shows  examples 
of  one-dimensional  block  pulse  functions. 
Two-dimensional  block  pulse  functions 
may  be  writ  ten  as  products  blo(K,x) 
blo(C,  y).  An  example  of  two-dimen- 
sional block  functions  is  shown  in 
figure  2. 

Block  pulse  functions  are  orthogonal, 
which  is  defined  as  (wo  functions  f(i,x)and 
Kj,x)  with  the  variable  x and  the  parame- 
ters i or  j in  the  interval  -Vi  < x < +'A  it  the 
integral  +'/2 

/ f(i,x)  f(j,x)  dx 

-Vi 

is  zero  for  i I i and  non-zero  for  i = j.  The 
functions  are  n irmaliz.ed  and  called  ortho- 
normal if  the  integral  equals  I for  i = j. 
from  figure  I,  the  integral  condition  for 
orthonomaly  is  verified  for  block  pulse 
functions. 
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Sets  of  orthogonal  functions  arc  suitable  for  decomposing  a given  function  in 
its  spectral  components  of  the  orthogonal  function.  The  most  commonly  known 
case  of  spectral  decompositions  is  the  lourier  analysis,  where  periodic  functions 
are  decomposed  into  an  infinite  serieso!  sine  and  cosine  functions  and  nonperiodic 
functions  are  decomposed  into  a lourier  integral  that  represents  the  continuous 
spectral  density  distribution  from  minus  to  plus  infinity.  Correspondingly,  the 
block  pulse  spectral  component  ol  a function  Itx)  existing  between  -14  and  + 14 
may  be  expressed  by 

+14 

Bin  = / f(x)  bio  (m,x)  dx 

-1/2 


If  an  image  of  the  dimension  a x b is  being  analyzed  by  a block  pulse  function, 
a is  scaled  to  lit  into  the  interval  -14  < x < 14  and  h into  -14  ( y ( 14.  The  two-dimen- 
sional spectral  component  B|nn  is  expressed  by 

■H/2 

B =J  / f(x,y)  blo(m,x)  b)o(n,y)  dx  dy 
mil  -Vi 


where  B)n  represents  the  average  transparency,  or  brightness,  of  the  image  in  tile- 
area  element  designated  by  the  block  pulse  value  bio  (m,x)  blo(n,y).  Ihe  spatial 
arrangement  and  characteristics  of  the  spectral  components  and  further  signal 
processing  may  be  used  to  identify  and  extract  topographic  terrain  features. 


Walsh  functions  are  another  example  of  a system  of  orthonormal  functions.  One- 
dimensional  Walsh  functions  are  expressed  by  wal  (k,x),  where  k represents  the 
order,  and  x the  variable.  These  functions  assume  the  two 
Functions  values  +1  and  - I only,  which  makes  them  comparatively  simple 

as  the  block  pulses  and  leads  to  simple  equipment,  flic  Walsh 
functions  are  defined  in  the  interval  -14  ( x < + 14  and  are  zero  everywhere  else. 
I hey  are  global  functions,  that  is,  the  functions  have  defined  values  over  the  entire 
interval.  I lie  shape  of  the  Walsh  function  is  more  complicated  than  the  block  pulse 
function  and  resembles  somewhat  the  sine  and  cosine  functions.  In  some  cases,  the 
Walsh  functions  are  periodic  within  the  interval;  in  other  cases,  they  are  nonperi- 
odic. The  functions  wal  (k,x)  may  be  defined  by  the  difference  equation 

wal  (2K+p,x)  = (-])IK/2I+P[wal  (K,2x+!4)  +(-l  )K+Pwal  (K,2x-'/2)| 


The  expression 


| K/2 1 means  the  largest  integer  less  than  or  equal  to  K/2. 
K =0,  I,  2,  3 n P = 0,  I 

1 1 for  -14  < x < 14  | 

and  wal  (o,x)  = 1 . . . | 

M)  lor  x < -14,  x > 14’ 
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I 


1 

I 

I 

' 
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In  the  following,  the  walsit  functions  wal  (l,x)  to  wal  (5,x)  are  developed  using 
the  generating  formula: 

K - 0 p = 1 

wal  ( 1 ,x)  - -wal  (0,  2 x + 14)  + wal  (0,  2x  - 14) 

K = 1 p = 0 

wal  (2,  x)  = wal  ( 1 , 2x  + Vi)  -wal  ( 1 , 2x  - 14) 

K = 1 p=  1 

wal  (3,  x)  ! -wal  (1 . 2x  + 14)  -wal  (1 , 2x  -14) 

K = 2 p = 0 

wal  (4,x)  = -wal  (2,  2x  + 14)  -wal  (2,  2x  - 14) 

K = 2 p = 1 

wal  (S.x)  = wal  (2,  2x  + 14)  -wal  (2,  2x  -14) 


figure  .3  shows  the  first  10  Walsh  functions  wal  (K,x).  All  of  these  functions  meet 
the  orthogonality  condition.  I here  fore,  any  continuous  function  defined  in  the 
interval  -Vi  ( x ( 'h  can  be  decomposed  into  spectral  components  of  a Walsh  function. 


I wo-dimensional  Walsh  functions  are  written  as  products  wal  (K.x)  (wal  fr.y) 
figure  4 shows  the  two-dimensional  Walsh  functions  from  wal  (0,x)  wal  (0,y) 
to  wal  (7.x.)  wal  (7,y.i.  I he  two-dimensional  spectral  Walsh  component  W^j, 
ol  the  function  f ( x,y ) is  +'/r 

W n = //  ft x,y ) wal  (m,x)  wal  (n,y)  dx  dy 

-Vi 


V*«  Of 


“1 

1 


**•  (2 
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if  1 < x ,y ) is  an  image  function  and  -'M  x V/i 
and  -'/A  y (Zi  is  scaled  to  match  the  dimen- 
sions of  the  image,  the  spectral  components 
W are  representative  of  the  entire  image 
rather  than  an  area  element  as  in  the  case  of 
spectral  block  pulse  components. 
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I inure  4.  I xjinplc  f >t  I wo-Dimcmional  Walsh 
I unctions. 
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Haar  functions  assume  the  values  of  +1,  0,  anil  -I.  The  first  few  functions  are 
shown  in  figure  5.  file  first  two  llaar  functions  are  identical  with  the  first  few 
Waish  functions.  The  one-dimensional  llaar  function  is  ex- 
pressed hy  liar  ( K,K,x ).  It  is  hart  0,  0,  x)  = wal  (0,  x ) 

liar  (0.  I . x ) = wal  (0.  x ) 

File  generating  process  for  Haar  functions  is  as  follows: 


llaar 

Functions 


I lie  first  two  Haar  functions  are  defined  as  the  first  two  Walsh  functions.  T he 
function  liar  < 0.  I,  \)  is  then  squeezed  into  half  the  interval  and  shifted  to  yield 
liar  (I.  1.x)  and  liar  (1,2,  x).  Kacli  one  of  the  new  functions  is  squeezed  again 
in  half  ol  its  space  and  shifted  to  yield  liar  (2,  1,  x),  liar  (2.  2,  x),  hat  (2,  3,  x), 
and  liar  (2,4.x).  I Ins  process,  if  continued,  generates  the  set  of  all  I laar  functions. 


fntirely  new  sets  of  discrete  orthogonal  functions  fn  can  he  generated  by  using 
Walsh  functions,  and  by  using  the  Haar  process  of  squeezing  to  one-half  and 
slid  ting.  I or  example. 


C Ct  f 


* v L 


't **• 


h*’ 


2*/ 


ha'  » v 1 


wal  (0,  x)  = f | (0,  0,  x) 

or  wal  ( D,  x)  = f-)  (0, 0,  x) 

wal  ( 1,  x)  = 1 ! (0,  1 . x) 

wal  (I.x)  1 A ( 0,  1 , x) 
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f|  (2,4.x) 

la  <2,  3,  X) 

— 

la  (2.  4,  x) 

^ a 

Two-  dimensional 

1 laar  and  1'  functions 

I ipurc  5.  I xatuplcs  ol  llaar  I unctions. 


are  expressed  as  products 

liar  ( m,  i,  x ) liar  ( n,  j,  y ) or 
l‘n  (in,  i,  x)  ln  (n,  j,  y)  respectively. 

As  witli  block  pulse  and  Walsh  functions, 
a spectral  decomposition  can  occur  with 
the  image  function  decomposing  into  a 
series  of  llaar  or  I components. 


A matrix  M is  called  an  oilhogonal  matrix  if  the  inverse  M * and  the  transposed 
M are  proportional  to  each  other.  A lladamard  matrix  is  an  orthogonal  matrix 
with  elements  t|  and  -I  only.  Il  is  convenient  to  write  instead 
<d  H and  -I  simply  + and  -.  Ihe  following  are  the  lladamard 
matrices  of  the  rank  I through  4: 


lludanuird 

Matrices 
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I 

I 


n,  = M 


'41 


+ + + + 
- - + + 
- + + - 
+ - + - 


II 


42  “ 


i+  + + - 
!+  + - + 

1+  - + + 
L + + + 


lladamard  matrices  of  higher  rank  can  he  generated  by  Kronecker  products,  litis 
process  is  explained  by  the  following  example: 


lh®ll 


+ + + + 

Ms  Ih 

+ - + - 

llr-lls 

+ + - - 

+ - - + 

because  the  rows  and  columns  of  the  lladamard  matrices  are  orthogonal  sets, 
they  can  be  used  to  decompose,  for  example,  a two-dimensional  image  function 
into  its  spectral  components.  The  orthogonal  condition  for  the  rows  of  an  llada- 
mard matrix  having  [sj-  elements  is  N 

V ;,.k  ;'jk  = No,,  Hv1,0'1  ' ' 

“ 11  'N  lor  i = j 

k = I 


I he  row  elements  ajk.K=  1,2...  .,N  may  be  considered  to  represent  a discrete 
function  that  assumes  the  values  I and  I within  the  interval  0<  x <N.  The  discrete 
function  may  be  expressed  by  a (i.x)  where 


a(  i.x ) = a,  j for  0 < x < I 
at  i.x ) = a,  •>  for  I < x ( 2 


I he  same  procedure  may  be  applied  for  the  column  elements  a,,  ^ i 1,2,3 N 

to  produce  a discrete  function  a(k.y). 

I he  corresponding  two-dimensional  spectral  components  for  ft  x,y)  are  expressed 
by 

N 

Itx.y)  at  i.x ) atk.y)  dxdy. 

() 


10 


I 


> 


i 

I 


I 


A shill  register  consists  of  n-consccutive  binary  storage  elements  called  stages  and 
of  a feedback  network  f igure  (i  shows  a shift  register  with  10  stages.  The  binary 
values  (e.g.  two  discrete  voltages)  that  can  be  stored  in  each 
stage  may  be  defined  as  value  /.ero  (0)  and  value  one  ( + 1). 


Discrete  Functions 
by  Shift  Registers 


H 


The  feedback  network  feeds  the  values  of  some  designated 
stages  into  a modulo  2 adder.  Depending  on  the  input  values  of  the  adder,  the 
output  of  the  adder  will  assume  the  value  0 or  +1.  The  output  of  the  adder  is 
applied  to  the  input  of  the  shift  register.  The  shift  register  is  operated  by  applying 

a clock  pulse  to  each 

-j  stage,  which  causes  the 

value  of  each  stage  to  be 
shifted  to  the  next  stage 
and  the  output  of  the 
feedback  network  to  be 
shifted  to  the  first, storage 
of  the  shift  register.  The 
output  of  the  shift  register 
7/J~r  is  then  a discrete  sequence 
of  0 and  1 . The  output 
series  may  consist  of  up 
to  2n  -1  bits  for  a n-stage 
register,  depending  on 
the  section  of  the  feed- 
back network.  The  word 
bit  refers  to  the  output  of 
the  shift  register  which  is 

generated  by  a one-clock  pulse.  Practically  all  discrete  functions  including  multi- 
value functions  may  be  generated  by  a suitable  set  of  shift  registers,  including 
switchable  feedback  connections  and  appropriate  logic  circuitry.  This  capability 
can  be  used  to  generate  many  reference  signals,  in  addition  to  orthogonal  se- 
quences, that  can  be  used  for  correlation  with  signals  derived  from  images. 


Figure  6.  10-Stage  Shift  Register  with  Feed  Hack  Network 


APPLYING  DISCRETE"  Discrete  function  technology  can  be  used  in  three 
FUNCTION  areas  of  the  topographic  sciences:  ( 1 ) in  aerial  photo- 

TECHNOLOGY  TO  graphs  of  topographic  terrain,  (2)  in  representing  the 
TOPOGRAPHIC  earth’s  gravity  field  by  Walsh  functions,  and  (3)  in  using 

SCIENCES  ‘I  iscrete  electro-magnetic  wave  forms  for  remote  sensing. 

Discrete  functions  can  be  electronically  generated  by  logic  switching  circuitry 
using  integrated  components.  I he  switching  speed  of  these  components  (T TL  and 
I CL)  is  typically  a few  nanoseconds  and  is  expected  to  increase  by  the  factor  100.' 
Therefore,  discrete  functions  may  be  generated  at  a tremendous  rate,  e.g.  up  to 
several  millions  per  second 

' Amin.  “High  Performance  iosephson  Circuits  Described  by  1HM  Researchers, "/FFF  Computer,  vol.lt 
No.  5,  p.  102,  May,  1978. 
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I he  area  elements,  or  pixels,  of  a black-and-white  photo  are  characterized  by  their 
pray  shades  or  degrees  of  transparency.  Hence,  the  image  can  he  described  as  a 
Filters  for  two-dimensional  gray  shade  function  expressed,  for  example, 

Two-Dimensional  in  ( artesian  or  polar  coordinates.  Topographic  and  terrain 
Image  Signals  features  recorded  by  a photographic  plate  or  film  are  recog- 

nizable by  their  shape,  contours,  and  gray  shade  distribution. 
Because  each  feature  has  its  characteristic  spectral  components  for  every  set  of 
orthogonal  functions,  spectral  components  of  features  may  be  useful  for  recog- 
nition, identification,  and  extraction  of  cartographic  features. 

Two-valued  discrete  functions  are  particularly  suited  for  development  of  spectral 
components  because,  they  can  easily  be  generated  with  off-the-shelf  components, 
furthermore,  electro-optical  devices  are  now  being  developed  that  can  be  used  for 
optical  generation  of  two-dimensional  discrete  functions.  The  one  value  0,  or  - I 
is  facilitated  by  dark  or  opaque  area  elements,  and  the  other  value  +1  by  light  or 
transparent  area  elements.  Four  types  of  electro-optical  filters  for  generating 
discrete  functions  are  described. 

Mematic  liquid  crystal  filters.  Liquid  crystals  are  organic  substances  whose  state 
1 of  matter  is  between  the  solid  crystaline  state  and  the  isotropic  liquid  state. 

The  liquid  crystals  are  classified  into  various  groups.  One  group,  for  electro-optical 
devices  is  the  nematic  liquid  crystals.  In  the  nematic  liquid  crystals,  the  molecules 
are  rod  shaped  and  arrange  themselves  parallel  in  volume  elements  containing 
typically  10^  to  I0<’  molecules.  The  volume  elements  have  characteristics  similar 
to  small  anisotropic  crystals.  For  electro-optical  applications,  many  of  the  volume 
elements  are  oriented  uniformly  by  putting  the  nematic  liquid  crystal  substance 
between  two  plane  glass  plates,  which  are  kept  at  a distance  of  10  to  100  micro- 
meters by  supports,  such  as  films,  threads,  etc.  The  inner  sides  of  the  glass  plates 
are  plated  with  transparent  electrodes  (Sn  0s,  In  0^)  that  make  it  possible  to 
effect  the  liquid  crystal  between  the  plates  by  electrical  fields.  I'lie  transparent 
electrodes  can  be  arranged  in  a pattern  of  parallel  stripes  isolated  by  very  small 
gaps  front  each  other.  To  make  an  optical  filter,  the  glass  plates  are  arranged  so 
that  the  stripes  of  the  one  plate  (x-direction)  are  perpendicular  to  the  stripes  of 
the  other  plate  (y-direction).  By  applying  voltages  on  crossed  stripes,  the  area  in 
between  can  be  made  either  transparent  or  opaque.  By  applying  the  values  of 
discrete  functions  electronically  to  the  x-stripes  and  y-stripes,  optical  filters 
representing  two-dimensional  orthogonal  discrete  functions  are  realized. 

Two  effects  of  nematic  liquid  crystals  can  be  exploited  for  the  realization  of 
optical  filters.  One  effect,  called  the  Dynamic  Scattering  (DS),  requires  non- 
polarized monocromatic  light  for  illumination.  The  other  effect,  called  De- 
formation Aufgerichteter  Phase  n (IMP)  requires  polarized  monochromatic  light 
lor  illumination.  The  contrast  ratio  betweeen  transparency  and  opaqueness  for 
the  DS  effect  is  1:50,  and  for  the  DAP  effect,  1:1000.  The  total  response  time, 
that  is  the  time  needed  for  switching  from  one  pattern  to  another,  requires 
approximately  3 microseconds  for  the  DS  effect  and  approximately  6 nano- 
seconds for  the  DAP  effect. 
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Photoferroelectric  (PFE)  ceramics.  Nonvolatile  images  can  he  stored  in  phase 
lead  lanthanum  zirconate  titanate  (I’LZT)  ceramic  by  using  near-ultraviolct(UV) 
light  and  switching  the  ferroelectric  polarization  through  a portion  of  the  hy- 
steresis loop.  The  images  are  stored  as  spatial  distribution  of  scattering  centers 
and  as  surface  deformations.  Erasure  of  stored  images  is  achieved  by  using  near- 
IJV  light  anil  switching  the  ferroelectric  polarization  back  to  its  initial  state.  In 
addition  to  the  image  storage  capability,  PIT  ceramics  may  be  used  for  devices 
that  generate  optical  filters  consisting  of  patterns  of  transparent  and  opaque 
areas.  These  filters  could  be  used  in  a similar  way  as  liquid  crystal  filters.  The 
advantage  of  PI  l.-PLZT  ceramics  is  that  they  are  solid  state  materials  and  can 
be  made  in  sheets  as  thick  as  0.5  millimeters.  However,  certain  problems  may 
arise  because  of  the  slow  total  response  time  and  possible  UV  absorption  in  the 
image  material. 

f igure  7 shows  an  optical 
configuration  for  spatial 
spectral  analysis.  The  par- 
allel light  beam  generated 
by  lens  1 is  modulated  by 
the  image.  A filter  function 
is  then  applied  by  the 
filter  to  the  modulated 
j light  beam.  Lens  2 focuses 

the  output  of  the  filter 
onto  a photo  cell.  The 
output  of  the  photo  cell 
is  a measure  of  the 
spectral  component  of 
the  image  in  *>*rms  of  the 
filter  function.  The  filters 
may  be  realized  by  liquid 
crystals  or  PFE  ceramics. 

Plasma  filters.  A new  type  of  display  technology  called  plasma  display  technology 
will  provide  a capability  to  generate  two-dimensional  light  distributions  that 
represent  discrete  functions.  Thin  sheets  of  glass  are  spaced  about  0.1  millimeter 
apart  to  form  a cavity  for  special  gaseous  mixtures.  Rows  of  x and  y conducting 
address  lines  are  deposited  on  this  glass.  The  x lines  are  on  the  one  side  of  the 
cavity;  the  y lines  are  on  the  other  side  and  perpendicular  to  the  x lines.  A thin 
dielectric  glass  film  is  placed  directly  over  those  lines,  and  a continuous  charge  of 
a.c.  voltage  is  applied  across  the  cavity  by  the  x and  y lines.  To  start  a discharge  in 
the  gas  at  any  point  where  an  x and  y line  intersect,  an  additional  voltage  pulse  is 
superimposed  that  is  large  enough  to  trigger  the  gas  discharge  at  this  point.  The 
gas  discharge  appears  as  a light  dot.  The  discharge  can  be  terminated  by  another 


/A  M til 


ligure  7.  Optical  Coni  juration  of  Spatial  Spectral  linage  Analyzer. 
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electrical  pulse  that  cancels  the  a.e.  voltage  at  the  line  intersection.  The  riots 
can  be  arranged  in  patterns  of  two-dimensional  discrete  functions,  such  as  Walsh 
functions.  This  light 
distribution  can  be  di- 
rectly applied  to  an  image 
to  produce  corresponding 

transforms  or  correlations.  \ > fttr.ro  , m 

The  transforms  can  be 
recorded  electrically.  The 
advantage  of  this  approach 
is  that  the  light  source 
and  optical  matrix  are 
combined  in  one  single 
device.  Such  a matrix 
may  be  called  an  active 
optical  filter,  figure  K 
shows  an  experimental 
setup  Using  a plasma 
filter  that  could  be  used 
for  pattern  recognition 
and  feature  extraction. 


>□- 


I i,  tin-  H I’riru  iplc  ( nnlifurutiun  <>l  .m  I Icctro-Opticul  Spatial  Spectrum  Analyzer. 


LEI)  arrays  as  filters.  I iglit  emitting  diodes  ( LI  U)  may  be  considered  as  building 
blocks  for  active  optical  lilters.  I hey  have  the  advantage  of  high  reliability  and 
very  short  switching  times.  The  si/.e  of  I I l)s  is  still  relatively  large  (about  I mini, 
and  the  problem  ol  manufacturing  large  two-dimensional  high  resolution  matrices 
has  not  yet  been  solved. 


Discrete  functions,  such  as  Walsh  functions,  can  also  be  used  to  construct  a binary 
sample  function  (BSI)  model  of  the  geopotential  and  to  recover  geoid  and 
geopotential  parameters  from  satellite  altimeter  data.  A 
comparison  of  the  performance  of  the  BSI  model  with  the 
performance  of  the  spherical  harmonics  sample  function 
( SI  IS  I- ) model  shows  that  the  BSI  model  is  belter  equipped 
to  recover  small  features,  typically  less  than  10  arc  degrees,  efficiently  and  accur- 
ately on  a side. 

A new  algorithm  for  rapid  spectral  analysis  of  satellite  altimeter  data  has 
been  developed  by  1<  I).  Brown2  that  reduces  the  computer  time  by  at  least  one- 
order  of  magnitude.  Ihe  algorithm  is  based  on  Walsh  functions  and  computes 
the  spectral  components  of  the  measured  data.  Ihe  algorithm  is  also  capable 


2Richard  I)  Drown.  “Cieopotcntial  Modeling  by  binary  Sample  functions,"  Catholic  University.  Washington, 
!).(.,  Doctoral  dissertation,  1974. 
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Representation  of 
the  Earth’s  (.ravity 
Field  by  Walsh 
Functions 
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<>l  handling  irregularly  spaced  data.  The  resulting  spectrum  retains  the  same  gen- 
eral shape,  amplitude,  and  location  of  peaks  as  the  spectrum  derived  from  the 
least-square  method.  If  N is  the  number  of  samples,  t|  the  required  processing 
time  for  the  Walsh  function  algorithm,  and  t->  the  processing  time  lor  the  least 
square  algorithm,  the  ratio  1 1 /t ->  is  kn  N/Nfcn  2.  l or  1 ,000  sample  points,  the  time 
t is  only  0.01  t->,  which  means  that  the  Walsh  process  is  100  times  faster  than  the 
least-square  process  for  this  sample  number. 


( urrently,  generators  can  be  built  that  can  produce  electromagnetic  discrete 
wave  forms  with  pulse  widths  as  small  as  I nanosecond.  However,  the  trend  is 
toward  shorter  pulse  widths  down  to  100  picoseconds  or  even 
to  10  picoseconds.  These  wave  forms  can  be  transmitted 
directly  by  an  array  of  radiators  that  serves  as  an  antenna 
system.  The  generator  and  antenna  systems  are  designed  into 
one  unit.  Because  the  current  that  is  fed  into  the  antenna 
elements  is  two  valued,  the  electrical  field  in  the  wave  zone  will  consist  of  positive 
anil  negative  Dirac  pulses. 


Applying  Discrete 
Function 
Technology  to 
Remote  Sensing 


Some  work  has  been  done  by  l-.T.  Ulaby3  at  the  University  ol  Kansas  in  which 
a sinusoidal  radar  spectrometer  was  focused  on  cropfields.  The  received  scattered 
energy  was  analyzed  as  a function  of  field  roughness,  soil  moisture  content, 
crop  type,  and  crop  height.  The  radar  used  a sinusoidal  carrier  of  4 to  18  (ill/ 
(gigahertz),  with  a I M (frequency  modulation)  of  approximately  400  Mil/ 
( megahertz).  I wo  of  Ulaby ’s  conclusions  are 


• 1 Although  soil  moisture  can  be  sensed  through  vegetation, 

the  sensitivity  of  radar  backscatter  to  soil  moisture  is  quite 
dependent  on  vegetation  characteristics  and  sensor  parameters. 

2.  Spectral  response  curves  indicate  that  lower  frequencies 
provide  more  information  on  soil  moisture  because  of  inher- 
ently better  penetrating  ability.  Angles  near  zennh  are  necessary 
to  accurately  estimate  soil  moisture. 

One  class  of  nonsinusoidal  waveforms  yields  more  information  than  sinusoids 
when  used  in  topographic  radar.  These  waveforms  are  sequences  of  short  pulses, 
particularly  positive  and  negative  Dirac  pulses,  where  an  amplitude  reversal  can  be 
easily  distinguished  from  any  phase  shift  even  in  a noisy  environment.  I his  effect 
also  exists  lor  sinusoidal  waves,  but  it  cannot  be  exploited  because  an  amplitude 
reversed  sinusoidal  function  looks  like  the  original  function  shifted  by  a hall 
period. 


*1  I I l|al>> , The  I tin  ts  nf  Soil  Moisture  ami  I’lant  Mor/’hology  on  the  Radio  Hat  kscatterer  from  li'ye  I a t loll . 
University  of  Kansas,  July,  1974. 
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Special  Walsh  functions  that  generate  two-valued  Dirac  pulses  have  been  used 
experimentally  lor  mine  detection,  detection  of  metallic  and  nonmctallic  pipes 
underground,  reflectors  in  noisy  environment,  etc.  Therefore,  ground  truth  data 
can  he  expected  to  be  more  accurately  obtained  by  using  discrete  function  radar 
rather  than  by  using  sinusoidal  pulsed  radar. 


FUTURE 
RESEARCH  AND 
DEVELOPMENT 


■ In  the  preceeding  section,  discrete  function  technology 
was  shown  to  have  a potential  for  applications  in  the 
topographic  sciences.  In  this  section, some  concepts  are 
presented  that  may  be  solutions  to  topographic  prob- 
lems that  are  now  significant  to  the  military. 
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figure  b shows  the  block 
diagram  ol  a concept  for 
image  analysis 
Cartographic  and  feature  ex- 
feature  traction  using 
Extraction  an  active  opti- 
and  cal  filter  (the 

Signal  filter  is  simul- 

Processing  taneously  the 
light  sourcefand 
a two-dimensional  sensing 
array  as  receiver.  I he  image 
is  placed  between  the  filter 
and  sensing  array.  Depen- 
ding on  the  resolution  of 
the  I liter  and  the  sensing 
array  spectral.  one  may 
obtain  decompositions  of 
the  entire  image  or  sections 


TWO  DlMENSIOAIAL 
■S£#3/Wa  At* AY 


1 /MAGE 

I igurc  9.  Block  iJugram  ol  A Concept  lor  image  Analyzer  and  I future  I xtractor. 


thereof  in  great  detail. 

Because  many  discrete  patterns  may  be  applied,  many  characteristic  spectra  may 
be  obtained  I he  required  function  generators  are  feasible  today.  However,  large 
active  filters  and  sensing  arrays  are  in  the  early  stages  of  development,  f urther- 
more. a detailed  analysis  ol  cartographic  features  is  required  for  suitable  two- 
dimensional  discrete  functions  that  produce  strong  responses  to  specific  features 
m terms  ol  spectral  lines  and  correlation.  I his  research  should  include  the  study 
ol  topographic  features  as  they  appear  under  various  sun  angles  for  clear,  partially 
cloudy,  and  cloudy  skies.  Also,  the  angle  of  the  object  relative  to  the  image  plane 
must  be  considered,  as  well  as  possible  shadows. 
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Sn  l;tr.  only  rectangular  matrices  have  been  considered  as  tillers.  One  idea  is 
to  include  other  curvilinear  coordinate  systems  as  matrices  to  express  discrete 
function  patterns.  figure  10  shows  two-dimensional  Walsh  lunetions  in  polar 
coordinates. 
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I i^iik'  Hi  I xampK  sol  Polar  Walsh  I unctions. 


Source  “Application  of 
Nonsinusoidal  Signals  to 
Cartographic  linage  Pro 
cessing,  ” Unpublished 
report  prepared  by  the 
Undersea  Research  < < »r p 
lor  the  U.S,  Army  Engi- 
neer I opographic  Labor- 
atories, Port  Helvoir , Vn 
ginia,  August  1 976. 


Another  idea  is  to  combine  the  photo  detector  lunetion  of  the  sensing  array 
with  tin-  Idler  lunetion  of  the  matrix.  I his  may  be  relatively  easy  to  accomplish 
loi  two  valued  Idler  lunetions.  One  value  would  be  represented  by  transparency, 
and  the  other  value  by  opaqueness  In  this  case,  the  opaqueness  may  be  realized 
by  disabling  the  corresponding  photo  cells  ol  the  array  utilizing  suitable  bias 
voltages.  I he  image  is  then  simply  projected  on  the  array.  Hits  also  gives  the 
possibility  ol  magnifying  sections  ol  the  image,  which  would  be  equivalent  to 
an  increased  resolution  ol  the  system. 


I he  success  ol  leal  lire  extraction  depends  to  a large  extent  upon  the  response 
time  id  the  components  lor  electro-optical  tillers  and  electronic  processing 
devices.  II  as  many  as  I to  100  million  signal  processing  procedures  can  be  ex- 
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cculcil  in  I second,  (lie  probability  lo  extract  topographic  feat  tires  may  be  high. 

I herclore,  the  Imal  liminatioii  is  presumably  the  speed  by  which  the  image  can  be 
moved  through  the  system  lor  image  analysis  and  leaturc  extraction. 

| he  Binary  Sampling  Model  is  a local  function  model  that  is  particularly  well 

suited  lor  puters.  I lie  strictly  local  nature  of  the  functions  also  allows  local 

geoid  solutions  lo  be  obtained  without  regard  to  the  lack  of  data  in 
(ifock-sy  areas  immediately  outside  the  local  solution  region.  I his  capability 
should  be  explored  when  the  geopotential  ol  limited  areas  ol  high 
military  importance  are  concerned.  Because  discrete  (unctions  may  also  be  cx- 
pressed  m terms  ol  surface  spherical  coordinates,  the  advantages  ol  expressing  the 
earth’s  gravity  field  by  these  functions  should  be  examined.  Figure  I I shows 
Walsh  (unctions  in  spherical  coordinates. 


Sensing 


Very  little  information  is  avail- 
able on  radar  type  instrumenta- 
„ . lion  using  discrete 

wave  forms.  I he  re- 
search performed  in 
this  area  is  primarily  con- 
cerned with  underwater  re- 
mote sensing  and  detection, 
usually  supported  by  the  Navy. 
However,  a remote  airborne 
sensing  radar  operating  with 
discrete  function  wave  forms 
would  provide  more  and  better 
vegetation  and  terrain  infor- 
mation than  sinusoidal  pulsed 
4 

radar.  Plans  are  being  made 
to  perform  experiments  with 
discrete-pulsed  radar  to  in- 
vestigate the  return  from  vege- 
tation and  terrain. 
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figure  1 1.  examples  of  Spherical  Walsh  functions. 
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■ On  the  basis  of  this  study,  the  following  conclusions  are 


1 . Discrete  function  technology  can  be  applied  to  at  least 
three  areas  of  the  topographic  sciences,  namely  to  image- 
analysis  for  cartographic  and  terrain  feature  extraction, 
to  geopotential  representation(and  to  remote  sensing. 

2.  Many  new  families  of  discrete  functions  can  be  derived  from 
the  Hadamard  matrices  by  means  of  the  Haar  process  and  by 
shift  registers.  These  functions  may  be  used  for  spectral  de- 
composition of  cartographic  features,  if  the  set  of  functions 
is  orthogonal,  or  for  correlation. 

3.  Currently,  hardware  components  are  available  to  verify 
concepts  of  experimentation,  which  includes  high  speed 
switching  circuitry,  nematic  liquid  crystals.  PM  -PLZT  ceramics, 
plasma  filters,  and  LI  I)  and  CCD  arrays. 

4.  High  speed  image  analyzers  and  feature  extractors  would 
require  component  response  times  of  approximately  1 nano- 
second to  complete  up  to  approximately  30  images  per  second. 

5.  Very  little  is  known  at  this  time  about  specific  spectral 
decompositions  of  topographic  features. 

<>.  In  situ  signal  processing  can  be  accomplished  by  using 
remote  sensing  array  cameras  combined  with  discrete  function 
processors  and  discrete  function  signal  transmissions. 

7.  Reduction  of  satellite  altimeter  data  can  be  accomplished 
very  efficiently  by  using  a binary  sampling  function  model. 

8.  Local  portions  of  the  geopotential  may  be  advantageously 
expressed  by  the  Walsh  functions. 

9.  Topographic  terrain  radar  will  provide  more  ground  truth 
information  when  operated  with  discrete  wave  forms  rather 
than  with  sinusoidal' wave  forms. 
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